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Abstract. We use the NASA GEOS-5 transport model with 
tagged tracers to investigate the contributions of different re- 
gional sources of CO and black carbon (BC) to their con- 
centrations in the Western Arctic (i.e., 50-90° N and 190— 
320° E) in spring and summer 2008. The model is evaluated 
by comparing the results with airborne measurements of CO 
and BC from the NASA Arctic Research of the Composi- 
tion of the Troposphere from Aircraft and Satellites (ARC- 
TAS) field campaigns to demonstrate the strengths and lim- 
itations of our simulations. We also examine the reliability 
of tagged CO tracers in characterizing air mass origins using 
the measured fossil fuel tracer of dichloromethane and the 
biomass burning tracer of acetonitrile. Our tagged CO simu- 
lations suggest that most of the enhanced CO concentrations 
(above background level from CH 4 production) observed 
during April originate from Asian anthropogenic emissions. 
Boreal biomass burning emissions and Asian anthropogenic 
emissions are of similar importance in July domain wise, al- 
though the biomass burning CO fraction is much larger in the 
area of the ARCTAS field experiments. The fraction of CO 
from Asian anthropogenic emissions is larger in spring than 


in summer. European sources make up no more than 10% 
of CO levels in the campaign domain during either period. 
Comparisons of CO concentrations along the flight tracks 
with regional averages from GEOS-5 show that the along- 
track measurements are representative of the concentrations 
within the large domain of the Western Arctic in April but 
not in July. 


1 Introduction 

Midlatitude pollutants, including short-lived species such as 

aerosols, are often transported to the Arctic (Barrie, 1986; 

Weber et al., 2003). A series of internationally coordinated 

Arctic field campaigns under the umbrella of Polar Study 

using Aircraft, Remote Sensing, Surface Measurements and 

Models, of Climate, Chemistry, Aerosols, and Transport 

(POLARCAT), e.g., Arctic Research of the Composition of 

the Troposphere from Aircraft and Satellites (ARCTAS) (Ja- 

cob et al., 2010), Aerosol, Radiation, and Cloud Processes 

affecting Arctic Climate (ARCPAC) (Brock et al., 2011), 
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POLARCAT-France (Quennehen et al., 2011), POLARCAT- 
GRACE (Pommier et al., 2010), and Airborne Extensive Re- 
gional Observations in SIBeria (YAK-AEROSIB) (Paris et 
al., 2009), were conducted in 2008 to better understand the 
impact of pollution on the Arctic atmospheric composition 
and climate. Here we use observations from the NASA ARC- 
TAS field campaigns including ARCTAS Phase-A, which 
was based in Fairbanks, Alaska, USA, in April and ARC- 
TAS Phase-B (ARCTAS-B), which was based at Cold Lake, 
Alberta, Canada, in July to study the source attribution of 
pollutants transported to the Western Arctic. 

Most of the pollution found in the Arctic air is from the 
midlatitudes, and a large fraction is attributable to anthro- 
pogenic emissions. It is necessary to quantify the sources of 
pollution and identify the transport pathways in order to de- 
velop effective control strategies. Previous studies identified 
European emissions as the main source of Arctic pollution 
(Barrie, 1986; Rahn, 1981; Raatz and Shaw, 1984; Quinn et 
al., 2007, 2008; Shaw 1995). Emissions from Europe have 
declined and anthropogenic activities in Asia have greatly in- 
creased in the past two decades (e.g., Streets et al., 2009). Re- 
cent investigations suggest that Asian pollution sources have 
become significant to Arctic pollution; their degree of im- 
portance is still being debated (Fisher et al., 2010; Koch and 
Hansen, 2005; Shindell et al., 2008; Singh et al., 2010; Stohl, 
2006). The uncertainty of the origin of Arctic pollution is 
partly due to complicated pathways for transport of air from 
midlatitudes to the Arctic (Fuelberg et al., 2010) and partly 
due to the uncertainties in estimating emission amounts and 
types in different regions. 

Based on the analysis of snow samples, Hegg et al. (2009) 
concluded that more than 90 % of the black carbon (BC) de- 
posited in the Arctic in spring is due to biomass burning, an- 
other important source of Arctic pollution. Fires from Russia, 
Eastern Europe, and Central Asia all have their imprints on 
the Arctic (Stohl et al., 2007; Warneke et al., 2009, 2010). 
However, unlike anthropogenic sources that are mostly from 
midlatitudes, local boreal forest fires occur in places much 
closer to or even within the Arctic Circle, making important 
contributions to the pollutant loading in the Arctic. 

A series of ARCTAS-related studies were published re- 
cently to address the source contribution to the Western Arc- 
tic, particularly by anthropogenic and biomass burning im- 
pacts. By using trajectories, Harrigan et al. (201 1) studied the 
mean transport characteristics in ARCTAS-A from anthro- 
pogenic emissions. A global chemistry and transport model 
GEOS-CHEM was used to interpret the ARCTAS-A mea- 
surements and to investigate the origin of carbon monox- 
ide (CO) (Fisher et al., 2010), ammonium sulfate aerosols 
(Fisher et al., 2011), and black carbon (BC) and organic 
aerosol (Wang et al., 2011) in the Arctic during spring sea- 
son. These studies indicated that Asian anthropogenic emis- 
sions and biomass burning are important sources for the 
Arctic spring pollution. Other studies used aircraft measure- 
ments to investigate pollution emissions, transport, atmo- 


spheric composition and chemistry, and spatial variations 
of pollutants (Cubison et al., 2011; Hecobian et al., 2011; 
Kondo et al., 201 la; Liang et al., 2011; Matsui et al., 201 la, 
b; McNaughton et al., 2011; Shinozuka, 2011; Shinozuka 
et al., 2011; Singh et al., 2010). To date there has been no 
modeling study exploring source attribution for the Western 
Arctic during both the ARCTAS-A and ARCTAS-B periods 
with an integrated analysis of CO and aerosols. Our work 
addresses this shortcoming. 

This study focuses on assessing the impact of long-range 
transport of midlatitude anthropogenic pollution and boreal 
forest fire emissions, both local and remote, on the West- 
ern Arctic atmospheric composition during both ARCTAS- 
A and ARCTAS-B. We use the NASA Goddard Earth Ob- 
serving System model version 5 (GEOS-5) simulations with 
tagged CO from different source regions and source types 
(Bian et al., 2010) to quantify the contributions of conti- 
nental sources and biomass burning to Arctic regional pollu- 
tion. Model-simulated concentrations are evaluated and con- 
strained by aircraft observations, and the tagged CO for air 
mass origin identification is compared with the observed an- 
thropogenic tracer dichloromethane (CH 2 CI 2 ) (Moore et al., 
2004) and the biomass burning tracer acetonitrile (CH 3 CN) 
(Gouw et al., 2003; Li et al., 2003; Singh et al., 2003). The 
model is then used for source attribution for the entire West- 
ern Arctic domain to put the ARCTAS measurements into a 
regional context. 

Carbon monoxide is emitted from both anthropogenic and 
biomass burning sources. It is not subject to dry and wet re- 
moval processes but can be oxidized by reaction with the hy- 
droxyl radical (OH). The lifetime of tropospheric CO is on 
the order of months (e.g., Pan et al., 1995), making it a good 
tracer for long-range transport. Black carbon, an important 
agent in Arctic climate change (Flanner et al., 2007), has 
similar sources as CO but a much shorter lifetime (about a 
week) because it can be removed by dry and wet depositions 
(Koch and Hansen, 2005). The relative abundance of CO and 
BC indicates the age of the air mass in the Arctic. 

We describe the NASA DC -8 measurements and NASA 
GEOS-5 model simulations during the ARCTAS period in 
Sect. 2. We evaluate the GEOS-5 model simulations using 
ARCTAS measurements in Sect. 3. In Sect. 4 we use GEOS- 
5 simulations to investigate sources, transport, and trans- 
formation of Arctic pollutants during the ARCTAS periods. 
Conclusions and discussions are given in Sect. 5. 


2 Aircraft observations and GEOS-5 model 

2.1 CO and BC concentrations from aircraft 
measurements 

The measurements of CO, BC, CH 2 CI 2 , and CH 3 CN 
from the NASA DC -8 aircraft (http://www-air.larc.nasa.gov/ 
cgi-bin/arcstat-c) are used to evaluate the model simulation. 
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There were 1 1 flights in March-April 2008 over Alaska and 7 
flights in July over Canada. CO was measured using tunable 
diode laser absorption spectroscopy (TDLAS), which pro- 
vides fast-response (~1 s) and high-precision (=b 1 parts per 
billion by volume) measurements (Diskin et ah, 2002; Sachse 
et ah, 1987). Fast-response aircraft data were averaged to 
1 min intervals along the flight tracks (http://www-air.larc. 
nasa.gov/cgi-bin/arcstat- c). Aerosol BC was measured us- 
ing a Single Particle Soot Photometer (SP2) based on laser- 
induced incandescence (LII) (Kondo et al., 2011b; Moteki 
and Kondo, 2007). The absolute accuracy of the BC measure- 
ments is estimated to be ~ 10% (Kondo et al., 2011a). The 
SP2 measures individual BC particles with a time resolution 
of 1-5 min depending on aerosol concentration. Whole air 
canister samples were used for analysis of dozens of volatile 
organic compounds, one of which was CH2CI2. Upon return 
to the laboratory, the canisters were assayed for CH2C12 by 
electron capture detection and mass spectrometry. The accu- 
racy of the measurement is +10 % and the measurement pre- 
cision was +10% (Colman et al., 2001). CH3CN was mea- 
sured by proton-transfer-reaction mass spectrometry (PTR- 
MS; Wisthaler et al., 2002). The precision of the 1 min aver- 
age data used for the analysis presented herein is ± 25 % at 
lOOppt. 

2.2 Model description 

GEOS-5 is a global Earth system model, containing com- 
ponents for atmospheric circulation and composition, ocean 
circulation and biogeochemistry, land surface processes, and 
data assimilation (Rienecker et al., 2008). The model has 
72 hybrid vertical sigma levels which are terrain following 
near the surface and transition to pressure levels above about 
100 hPa with a model top at 0.01 hPa (about 85 km). The 
model was run at 0.5° x0.625° latitude x longitude horizon- 
tal resolution. The model is “replayed” from the Modern- 
Era Retrospective Analysis for Research and Applications 
(MERRA) meteorological analyses at the same spatial res- 
olution produced by the NASA Global Modeling and Assim- 
ilation Office (Rienecker et al., 2011). Every 6h the model 
dynamical state (winds, pressure, temperature, and humidity) 
is set to the balanced state provided by MERRA and then a 
6 h forecast is performed until the next analysis is available. 

A version of the Goddard Chemistry, Aerosol, Radiation, 
and Transport model (GOCART; Chin et al., 2002; Mian 
Chin et al., 2009) was implemented in GEOS-5 and used 
to simulate processes of sources, sinks, transport, and trans- 
formation for CO, SO2, and aerosols within the GEOS-5 
system. A description of the aerosol simulation in a previ- 
ous version of the GEOS modeling system (GEOS-4) was 
provided in Colarco et al. (2010). GOCART aerosols in- 
cluded dust, sea salt, sulfate, black carbon, and organic mat- 
ter, mixed externally (non-interacting). In this study GO- 
CART CO was decomposed into five tagged CO tracers, 
each designed to track an important emission type or loca- 


Table 1 . Lifetime (days) of BC against dry, wet, and total deposi- 
tions and the contribution of wet deposition to the total deposition. 



Lifetime (days) against 

Wet dep. 
contribution (%) 

April 

Dry dep. Wet dep. Total dep. 

24 7.2 5.5 

77% 

July 

29 7.7 6.1 

79% 


tion (Bian et al., 2010; Fisher et al., 2010; Fuelberg et al., 
2010; Stohl, 2006). The regions of the five CO source types 
are shown in Fig. la. The five CO sources include three an- 
thropogenic sources from North America (NAFF or FI, yel- 
low), Asia (ASFF or F2, orange), and Europe (EUFF or F3, 
green) and two biomass burning emission sources from non- 
boreal regions (NBBB or B 1 , blue-shaded area) and boreal 
regions (BOBB or B2, non-blue-shaded area). Figure lb-c 
shows the DC-8 flight tracks for all flights during ARCTAS- 
A and ARCTAS-B, respectively. 

Two chemical reactions are included in the CO simula- 
tion: chemical loss of CO via reaction with OH, and chem- 
ical production of CO from CH4 oxidation by OH. CO pro- 
duced from non-methane hydrocarbon oxidation is parame- 
terized as a direct emission (Bian et al., 2007). The OH field 
is prescribed using the results from a global chemistry and 
transport model simulation, the Global Modeling Initiative 
(GMI). The CH4 field is prescribed also from GMI based 
on the measurements at the worldwide NO A A Global Moni- 
toring Division (GMD) sites and distributed as a function of 
latitude (Bian et al., 2007). 

Black carbon is emitted as 80% hydrophobic and 20% 
hydrophilic (Liousse et al., 1996; Chin et al., 2002). Dur- 
ing the “aging” process, the hydrophobic BC is converted to 
the hydrophilic phase with an e-folding time of 1.1 to 2.5 
days (Cooke et al., 1999; Wang et al., 2011; Maria et al., 
2004). The hydrophilic BC is subject to wet scavenging that 
accounts for the scavenging in convective updrafts and rain- 
out/washout in large-scale precipitation (Giorgi and Camei- 
des, 1986; Balkanski et al., 1993). The wet scavenging coef- 
ficient of hydrophilic BC is set to be 0.4. Dry deposition of 
BC is calculated by a resistance-in-series scheme (Wesely, 
1989). 

We computed the GEOS-5 global lifetimes of BC against 
dry and wet deposition for April and July 2008 (Table 1). The 
lifetime of BC is 5.5 days in April and 6. 1 days in July, which 
is comparable to the annual mean value from the GEOS- 
Chem model for the same year (5.9 days, Wang et al., 2010). 
Wet deposition contributes 77 % and 79 % to total removal of 
BC in April and July, respectively - also close to the annual 
average value from GEOS-Chem (77 %, Wang et al., 2010). 
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Table 2. Global and regional emissions of CO and BC in April and July 2008. 



CO (TgCO mon - 1 ) 

BC (TgC mon -1 ) 


April 

July 

April 

July 

Anthropogenic 

57.7 a (51/59) b 

59.6 

0.45 

0.47 

NAFF (FI) 

7.3 

7.6 

0.032 

0.033 

ASFF (F2) 

28.8 

29.7 

0.25 

0.26 

EUFF (F3) 

5.4 

5.6 

0.046 

0.047 

Secondary production from NMHC oxidation 0 

9.7(84/9.4) 

9.9 



Open fires 

43.6(104/50) 

46.1 

0.63 

0.48 

NBBB (Bl) 

19.6 

33.5 

0.30 

0.42 

BOBB (B2) 

19.7 

8.0 

0.33 

0.060 

Secondary production from NMHC oxidation^ 

4.3 (10/4.9) 

4.6 



Biogenic 

21.9 (29/29) 

36.3 



Methane 

67 (71/71) 

75 



Ship 



0.0049 

0.0050 

Total 

190 (255/209) 

217 

1.1 

0.95 


a All numbers outside parenthesis are from this work. 

b Numbers inside parenthesis are from the GEOS-Chem prior (first number) and optimized (second) emissions from Fisher 
et al. (2010). 

c The emission is included in the anthropogenic emission total. 
d The emission is included in open fire emission. 

2.3 Emissions of CO and BC 

The emissions for the CO and aerosol simulations are as de- 
scribed in Bian et al. (2007) and Colarco et al. (2010), ex- 
cept as updated here. Anthropogenic emissions of CO are 
adopted from the Emission Database for Global Atmospheric 
Research (EDGAR) 2000, overwritten regionally by Envi- 
ronmental Protection Agency (EPA) 1999 National Emis- 
sion Inventory over USA, European Monitoring and Eval- 
uation Programme (EMEP) over Europe, Big Bend Regional 
Aerosol and Visibility Observational (BRAVO) over north- 
ern Mexico, Streets 2000 over Southeast Asia, and Streets 
2001 over China (Yuan et al., 2012). Anthropogenic emis- 
sions of BC (including ship and aircraft emissions) are from 
the compilation of A2-MAP emission that is available for 
global model simulations (Diehl et al., 2012), such as models 
in the international project of Aerosol Comparisons between 
Observations and Models (AeroCom). The biogenic source 
of CO is prescribed based on results from the Global Mod- 
eling Initiative (GMI) global chemistry and transport model, 
which uses the Model of Emissions of Gases and Aerosols 
from Nature (MEGAN, Guenther et al., 2006) for biogenic 
emissions of non-methane hydrocarbon (NMHC) as a CO 
source. CO from methane oxidation is included as well (de- 
tails in Bian et al., 2007). 

Biomass burning emissions of CO and BC are from the 
Quick Fire Emission Dataset (QFED) version 2.4. The QFED 
emissions are based on the fire radiative power (top-down) 
approach to derive the amount of combusted biomass (i.e., 
dry mass burned or DM) by calibrating initially against 
global monthly mean emissions of Global Fire Emission 
Database (GFED) version 2 (Darmenov and da Silva, 2013; 


Petrenko et al., 2012). The fire radiative power (FRP) and lo- 
cation of the fires are obtained from the Moderate Resolution 
Imaging Spectroradiometer (MODIS) Level 2 fire products 
(MOD 14 and MYD14) and the MODIS Geolocation prod- 
ucts (MOD03 and MYD03). Species-specific emission fac- 
tors are then applied to DM estimates to obtain tracer gas 
and aerosol emissions. For aerosol species, the QFED emis- 
sions are further calibrated by employing biome-dependent 
enhancing factors based on the constraint on model AOT 
with MODIS AOT for four biomes under our current as- 
sumptions of particle properties (size and refractive index, 
and hence mass extinction efficiency, MEE). For example, 
the enhancement factors of BC are ranging from 1.8 (tropi- 
cal biome) to 4.5 (extra-tropical biome). 

Table 2 summarizes the emissions of CO and BC over 
the globe and in the five tagged region domains. Most of 
the CO is emitted from anthropogenic and biomass burn- 
ing sources, with anthropogenic emissions being about 30 % 
higher than biomass burning. Regionally, the important emis- 
sions are anthropogenic emission from Asia (ASFF) and 
biomass burning emissions from mid to high latitudes in 
the Northern Hemisphere (NH) (BOBB) in both months. 
Further analysis indicates that high BOBB occurs primar- 
ily over Russia/Kazakhstan in April and over North Amer- 
ica in July. CO emission in July is higher than in April by 
15 % in total. Fisher et al. (2010) studied CO over the Arctic 
during ARCTAS-A using the GEOS-CHEM model. When 
compared to ARCTAS-A observations, they found that their 
sources needed to be adjusted to optimize the agreement with 
the observations. The GEOS-5 CO emissions are comparable 
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Fig. 1. (a) Regions for tagged CO emissions. In this work we use 
3 tagged fossil fuel CO tracers and 2 tagged biomass burning trac- 
ers: FI. North American fossil fuel (NAFF, yellow), F2. Southern 
Asian fossil fuel (ASFF, red), F3. Europe fossil fuel (EUFF, green), 
Bl. Non-boreal biomass burning (NBBB, light-blue shaded area), 
and B2. Boreal biomass burning (BOBB, non-light-blue area), (b) 
and (c) NASA DC-8 flight tracks during ARCTAS-A (April, based 
at Fairbank, USA) and ARCTAS-B (July, based at Cold Lake, 
Canada), respectively. 


with the optimized GEOS-Chem CO emissions in April (see 
GEOS-Chem emissions in Table 2). 

BC is also emitted from anthropogenic and biomass burn- 
ing sources. In April the biomass burning emission is 30- 
40 % higher than anthropogenic emission, but in July these 
two sources are roughly comparable. The most important 


sources of BC are similar to those of CO, i.e., ASFF and 
BOBB, in both April and July. 

Uncertainty in the emissions of CO and BC is quite large, 
particularly in specific regions. Shindell et al. (2008) inves- 
tigated pollution transported to the Arctic using 17 global 
models. The difference of the emissions of CO and BC 
among the models are as great as a factor of two between 
the minimum and maximum over the four polluted regions 
of East Asia, South Asia, North America, and Europe in their 
baseline simulation. In comparison, emissions of CO and BC 
over East and South Asia in GEOS-5 are close to the maxi- 
mum value, while the emissions over North America and Eu- 
rope are close to or even slightly below the minimum value of 
the multi-model range in Shindell et al. (2008). Considering 
that their baseline simulation was for the year 2001 and the 
emission trends differ over the regions (i.e., increase ^20 % 
of anthropogenic emission over Asia and decrease ^40% 
over North America and Europe from 2001 to 2008), GEOS- 
5 emissions should be consistent with the other models. 


3 Model evaluation and analysis of ARCTAS data 
3.1 CO and BC concentrations 

Vertical distributions of CO volume mixing ratios from DC-8 
measurements and GEOS-5 simulation are shown in Fig. 2a 
and b for ARCTAS-A and -B in April and July 2008, respec- 
tively. The model output was sampled at the time, latitude, 
longitude, and altitude of the 1 min merged DC-8 measure- 
ments, with both the model and data values further averaged 
in 1 km vertical bins. The mean and standard deviation of 
the data are shown as red bars in Fig. 2. The modeled CO 
from other than the aforementioned five tagged categories is 
designated as “Other CO”. This “Other CO” comes mainly 
from the CO produced by CH4 oxidation, which typically 
accounts for about 30 % of global CO source (e.g., Arellano 
et al., 2006; Bian et al., 2007; Duncan et al., 2007). 

The April campaign (ARCTAS-A) was conducted in 
Alaska, focusing on the long-range transport of midlatitude 
pollutants. This long-range transport can be seen by a large 
fraction of CO from ASFF in Fig. 2a. Together, the midlati- 
tude anthropogenic emissions from Asia, Europe, and North 
America account for nearly 50% of CO along the DC-8 
flight track in April, with the ASFF contributing most (about 
25 %). 

The July campaign (ARCTAS-B) was based at Cold Lake, 
Canada. The areas covered by the July campaign experienced 
strong influence from local fires so that model calculated CO 
from BOBB is extremely high (> 60 % of total five CO tags) 
from the surface up to 3 km (Fig. 2b). CO transported from 
Asia peaks at high altitudes (~9 km) contributing to about 
half of the CO there. In both campaigns the background 
CO, mostly produced from CH4 oxidation, is high (about 
1/3 of total CO as shown by the “Other CO”). The general 
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Fig. 2. Vertical distribution of the CO volume mixing ratio (ppbv) 
from DC-8 measurements and the GEOS-5 simulation when the 
GEOS-5 model results are sampled by all flights for April (a) and 
July (b). The DC- 8 CO from along the flight tracks is shown by 
the solid thick red line with standard deviation shown by horizontal 
bars. GEOS-5 CO is shown by the thick black line for total and by 
the color-shaded areas for the corresponding tagged components. 
The “Other CO” refers to the global CO other than the five tagged 
COs defined in Fig. la and it is mostly attributed to the CO chemical 
formation from methane oxidation. 


April July 




Fig. 3. Vertical distribution of BC mass concentration (ngm -3 ) 
from DC-8 measurement and GEOS-5 simulations when the GEOS- 
5 model results are sampled by all flights for April (Fig. 2a) and 
July (Fig. 2b). Four model sensitivity experiments were conducted 
to test the different biomass burning emissions (green line for QFED 
2.4 and red line for Modified QFED 2.4) and BC aging e-folding 
timescales (solid line for e-folding time 2.5 days and dash line for 
1.25 days). The thick lines (i.e., green solid line in April and red 
dashed line in July) are the “best” simulations used in the discus- 
sion of source attribution. Correlations between the measurement 
and the best model simulation are given. 


consistency between model and observations on CO concen- 
trations and vertical profiles suggests that the GEOS-5 is ca- 
pable of reproducing the average levels of observed CO. 

The model-measurement comparison for BC is shown in 
Fig. 3 with similar data binning as in Fig. 2. Four model 
sensitivity simulations are shown in the figure to test the 
range of contributions from biomass burning emissions and 
BC aging times (conversion time from hydrophobic to hy- 
drophilic). The green lines represent BC biomass burning 


emission from QFED v2.4, while red lines stand for modi- 
fied QFED v2.4 BC emissions where the enhancement factor 
over extra-tropical areas was set to 1 . The different line styles 
are associated with BC e-folding aging times, with solid lines 
for 2.5 days and dashed lines for 1.25 days. 

The changes of BC biomass burning emission and e- 
folding aging times both have significant impact on modeled 
BC vertical profiles in April. The best agreement comes from 
the simulation driven by standard QFED v2.4 emission with 
an e-folding aging time of 2.5 days (green solid line). On the 
other hand, in July the BC concentrations at lower altitudes 
(e.g., below 4 km) are predominantly controlled by the local 
biomass burning emissions as the DC -8 often flew through 
the biomass burning plumes. Figure 3b clearly suggests that 
the biomass burning emission from the standard QFED is too 
high (by a factor of 7). By removing the “enhancement fac- 
tor” from the standard QFED, the agreement between model 
and observation is much improved, even though the model is 
still 50-70% higher than the observations below 4 km. Al- 
though the change of BC aging time makes little difference 
on the average BC vertical profile in Fig. 3b, using a faster 
e-folding time allows the model to better represent the “back- 
ground” BC (i.e., BC concentrations below 10 ng m -3 ) when 
compared to measurements from each individual flight (see 
Fig. S4 in the Supplement). Therefore, using the ARCTAS 
data as guidance, we select the optimal parameters within 
the range of BC e-folding time (Sect. 2.2) and biomass burn- 
ing emission (Sects. 2.3 and 3.1) in model simulation: the 
adjusted (reduced) biomass burning emission in July, and a 
season-dependent BC aging time (slower in April than July) 
that possibly reflect the influences of humidity and tempera- 
ture on the BC aging process. 

More evaluations of model simulation are presented in the 
Supplement. CO and BC concentrations from model and ob- 
servation along each flight track during April and July are 
shown in Figs. S1-S4. A composite of the DC- 8 -measured 
CO mixing ratios compared to the total CO mixing ratio sim- 
ulated by the GEOS-5 model along the flight tracks is pre- 
sented in Fig. S5, along with a similar comparison for BC in 
Fig. S 6 . 

3.2 Air mass origin 

GEOS-5 uses 5 tagged CO tracers (Fig. la) to track contri- 
butions of different sources to the Western Arctic. To eval- 
uate the reliability of this model, we analyze the correlation 
of the modeled anthropogenic CO (i.e., FF, the total of ASFF, 
NAFF, and EUFF) with observed dichloromethane (CH 2 CI 2 ) 
(Fig. 4a and b), a well-known fossil fuel tracer (Moore et al., 
2004), and the correlation of the modeled biomass burning 
CO (i.e., BB, the total of BOBB and NBBB) with observed 
acetonitrile (CH 3 CN) (Fig. 4c and d), a well-known biomass 
burning tracer (Gouw et al., 2003; Fi et al., 2003; Singh et al., 
2003). Both tracers have a lifetime on the order of 5 months. 
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Fig. 4. Top: scatter plots of observed CH2CI2 (pptv) and modeled fossil fuel tagged CO (ppbv) from Asia, North America, and Europe in 
April (a) and July (b). Bottom: scatter plots of observed CH3CN (pptv) and modeled biomass burning tagged CO (ppbv) from boreal and 
non-boreal regions in April (c) and July (d). The results are colored according to the measurement altitudes. Each point represents median 
values of model results and observations within 10 min. 


The modeled FF CO is correlated with observed CH2CI2 
in both campaigns (i.e., both R > 0 . 6 ), as shown in Fig. 4 a 
and 4 b. This is direct evidence that the model can capture 
the air mass that is transported from the midlatitudes well. 
The linear fitting analysis indicates that the ratio of tagged 
FF CO and observed CH2CI2 is higher in July than in April. 
The CO chemistry does not support this since more CO is 
lost to its reaction with OH in July. To answer this a detailed 
analysis of the emission ratio between CO and CH2CI2 in 
every anthropogenic emission sectors is needed, but the study 
is beyond the scope of this paper. 

In contrast to the correlation between tagged FF and ob- 
served CH2CI2, the correlation between the tagged BB CO 
and the observed CH3CN is much more complicated (Fig. 4 c 
and d). Even though CH3CN is a widely used biomass burn- 
ing tracer, our analysis indicates that the CH3CN and BB CO 
do not always covary. There are clearly three branches in the 
CO-CH3CN relationship. Branch 1 shows covarying CO and 
CH3CN concentrations that predominantly result from their 
common biomass burning sources. Branch 2 shows an ele- 
vated CO but low CH3CN when sampling from Arctic ma- 
rine boundary layer (AMBL) air since CH3CN is taken up 
by the ocean. Branch 3 shows an elevated CH3CN but low 
CO in stratospheric air since CH3CN has a longer lifetime 
than CO. These behaviors are nicely confirmed by the alti- 


tude color coding. Branch 1 clearly controls the CO-CH3CN 
relationships in July as the DC-8 frequently sampled biomass 
burning air masses during ARCTAS-B. This 3 -branch behav- 
ior is shown better in the scatter plot of total CO and observed 
CH3CN, in Fig. 5 , where the total CO shown originates from 
both the observations and the model results. The characteris- 
tics of these different CO-CH3CN relationships are captured 
by the model (Fig. 5 b and d). 

To demonstrate the utility of air mass identification us- 
ing the tagged tracers or observed CH2CI2 and CH3CN as 
tracers, we show in Fig. 6 the flight on 8 July from Cold 
Lake, Canada, to Thule, Greenland (ARCTAS-B, flight 21 ), 
as a detailed case study. This flight encountered two high 
CO plumes: one over the US-Canada border (waypoints 1 - 3 ) 
and the other over the Arctic near Greenland (waypoints 8- 
10 ) (see Fig. lc for the flight route, the waypoints are evenly 
distributed along flight track). The simulated tagged CO con- 
centrations along the flight tracks are shown in Fig. 6a and 
the observed CH2CI2 and CH3CN in Fig. 6b. To better un- 
derstand the contributions from different source regions, cur- 
tain plots of the tagged CO along the flight track are shown 
in Fig. 6c-d, respectively, for its major sources BOBB and 
ASFF. 

Analysis of the atmospheric circulation revealed that the 
jet stream consisted of two segments over Asia. The southern 
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Fig. 5. Scatter plots of observed CH3CN (pptv) and DC-8 observed CO (ppbv) (a and c) and of observed CH3CN (pptv) and GEOS-5 
simulated total CO (b and d) in ARCTAS-A and -B campaigns. The results are colored according to the measurement altitudes. Each point 
represents median values of model results and observations within 10 min. 


segment passed over Central Asia and was near the clima- 
tological position, while the segment over extreme northern 
Asia was ~10ms -1 stronger than its weaker climatologi- 
cal counterpart (Fig. 16 in Fuelberg et al., 2010). Ten day 
back trajectories from the aircraft tracks show that both the 
southern and northern plumes came from Asia and its sur- 
rounding areas. More specifically, the southern plume orig- 
inated from Asia and Eurasia and the northern plume from 
Northeast Asia and the Labrador Sea (Fuelberg et al., 2010, 
http://fuelberg.met.fsu.edu/research/arctas/traj/traj/html). 

The tagged CO simulations, indicated by the colored lines 
in Fig. 6 a, show that the sampled CO concentration in the 
southern plume had comparable Asian industrial and boreal 
biomass burning pollution contributions, while the sampled 
CO in the northern plume was mainly due to Asian industrial 
pollution, consistent with the origin of the plumes analyzed 
by Fuelberg et al. (2010). The contribution and distribution of 
BOBB and ASFF can be identified more easily from Fig. 6 c- 
d. Although both emission sources contributed to the two 
plumes, their peaks are slightly displaced vertically. In par- 
ticular, for the northern CO plume, the Asian anthropogenic 
influence peaks at 300-400 hPa, higher than the 400-500 hPa 
peak of the calibrating initially against plume, so the flight 
captured mostly ASFF at its flight altitude. This layering of 
plumes, i.e., the ASFF plume just above the BOBB plume, 


was also observed by previous studies (Matsui et al., 2011a; 
Singh et al., 2010). 

In general, the observed CH 2 CI 2 and CH 3 CN have simi- 
lar variations along the flight track as the model ASFF and 
BOBB, respectively. This comparison further indicates that 
the tagged CO tracers adequately represent the origin of air 
masses. 

3.3 Age and origin of air masses diagnosed by the 
BC/CO ratio 

As noted earlier, CO and aerosol BC have common sources 
from anthropogenic and biomass burning but different re- 
moval and chemistry processes. Given these similarities and 
differences, it is possible to use correlations of BC and CO 
to estimate air mass ages and origins. 

The ratios of BC to CO from DC -8 measurements and the 
GEOS-5 model in the two ARCTAS campaign periods are 
shown in Fig. 7. The older the air is, the smaller the BC/CO 
ratio will be, simply because the lifetime of BC is much 
shorter than CO since BC is removed much faster through 
dry and wet scavenging than CO. The ratios in the figure are 
calculated from BC and total CO, but are sorted according to 
the dominant sources of tagged CO. Here the dominant CO 
refers to the CO tag having the greatest value, as well as its 
value being not less than 40 % of the total CO. We show only 
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Fig. 6. Top left: the five tagged CO mixing ratios during the ARC- 
TAS flight 21 on 8 July from Cold Lake, Canada, terminating in 
Thule, Greenland. Top right: the corresponding measured CH3CN 
and CH2CI2 during the same flight. Bottom left and bottom right 
are curtain plots for the tagged CO species BOBB and ASFF, re- 
spectively, with the white line representing the flight track. 
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Fig. 7. The statistical ratios between the mixing ratios of BC 
[ngsm -3 ] to CO [ppbv] over ARCTAS-A (top) and ARCTAS-B 
(bottom) campaigns. The blue box represents model ratio and yel- 
low box represents observed ratio. The ratios are sorted based on 
the dominant model CO tag in the measured pollution. 


Table 3. The ratio of BC / CO (Gg Tg 1 ) from emissions in differ- 
ent source regions and from Arctic air samples in July. 


BC(Gg)/CO(Tg) 

ASFF 

NAFF 

BOBB 

NBBB 

Emission 

14 

4 

8 

13 

Air sample (model) 

0.2 

2 

2 

2 

Air sample (measurement) 

0.1 

2 

1 

2 


Asia and other tropical regions such that CO and BC from 
NBBB had traveled a long distance to reach Alaska. In 
ARCTAS-B, the BC/CO ratio is the lowest in the ASFF air 
mass (~0.15 BC(ng sm _3 )/CO(ppbv)), but similar in NBBB 
and NAFF samples (~1.7 BC(ng sm _3 )/CO(ppbv)), reflect- 
ing the longer traveling distance of Asian air mass to Canada 
than the North American air mass that contains pollutants 
from NAFF and NBBB (Californian fire). Although the mean 
BC/CO ratio of NOBB from GEOS-5 is higher than that 
from the observation, model and measurement data are con- 
sistently showing that the BC / CO ratio has the largest vari- 
ation in BOBB (indicated by the highest standard deviation) 
because of the two distinct origins: one was from local Cana- 
dian boreal fire (fresh) and the other was transported from 
Eurasia region (aged). The BC / CO ratio in ASFF was larger 
in April than in July, which is at least partially due to more 
BC wet scavenging occurring during the Asia-Arctic trans- 
port in July (Matsui, et al., 201 la). 

We also do a first-order check for the ratio of BC / CO 
in ARCTAS-B by comparing their values with the ratios 
of BC / CO from emissions data. Theoretically, the former 
should be smaller than the latter since BC is lost faster than 
CO during transport. In July the emission ratios of BC / CO 
(Table 3) are in the 4-14 (Gg Tg -1 ) range for the five tagged 
CO tracers, which are larger than the BC / CO ratio in the air 
samples that are within 0.2-2 (Gg Tg - 1 ) in the tagged CO ap- 
proach and 0.1-2 (Gg Tg -1 ) from observations. The BC/CO 
ratio of the tagged ASFF air mass has the largest change from 
its emission to the target air sample. This is understandable, 
given that the emission from the ASFF is more remote and 
the large precipitation was reported along trajectories from 
Asian anthropogenic air to the Cold Lake, Canada, in July 
(Matsui, et al., 201 la). 


for those dominant CO tracers when there are more than ten 
1 min sample numbers having both BC and CO values, so 
there is no ratio calculated for cases of EUFF in both cam- 
paigns and NAFF in April. 

In ARCTAS-A the BC / CO ratio is the smallest for NBBB 
(~0.15 BC(ng sm -3 )/ CO(ppbv)), indicating that air coming 
from NBBB was the oldest (unit conversion from volume 
mixing ratio to mass concentration using standard temper- 
ature and pressure condition). During April the major NBBB 
emissions were from the low-latitudinal areas of Southeast 


4 Regional source attribution of CO for the 
Western Arctic 

Since CO and aerosols are not uniformly distributed and the 
aircraft measurements can only cover a small area in a lim- 
ited time frame, the question often arises: how representative 
are the findings from the ARCTAS field experiments for the 
larger Arctic region? Here we use the GEOS-5 model to es- 
timate the pollution sources in the Western Arctic domain to 
put the ARCTAS results into a regional context. 
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Fig. 8. Vertical distribution of the CO volume mixing ratio (ppbv) 
from GEOS-5 simulation when the model results are averaged 
over the Arctic region 50N-90N and 190E-320E for April (a) and 
July (b). The GEOS-5 CO is shown by the thick black line for total 
and by the color shaded areas for the corresponding tagged compo- 
nents. The “Other CO” refers to the global CO other than the five 
tagged CO species defined in Fig. la and it is mostly attributed to 
the CO chemical formation from methane oxidation. The CO mix- 
ing ratio from along the DC-8 flight tracks is also shown by the 
dotted red line with the standard deviation shown by horizontal bars 
for a reference. 


The vertical profiles of the CO volume mixing ratio av- 
eraged over the entire Western Arctic domain (between 50- 
90° N and 170-40° W) from GEOS-5 are shown in Fig. 8a 
(April 2008) and Fig. 8b (July 2008). The CO results from 
DC-8 measurements are plotted in dotted lines. The similar- 
ity between Figs. 2a and 8a in terms of CO concentrations 
and source attribution implies that the pollutants in April 
2008 were relatively well mixed and homogenously dis- 
tributed in the Western Arctic and the ARCTAS-A data pro- 
vide an adequate representation of the region, even though 
the flight planning during the mission often directed aircraft 
to capture the transported pollution plumes. In contrast, in 
July the fire plumes from BOBB in Canada were specifi- 
cally targeted during ARCTAS-B, which were highly inho- 
mogeneous in both space and time. The large difference be- 
tween Fig. 2b (flight track average) and Fig. 8b (regional av- 
erage) in total CO and source apportionment from GEOS-5 
clearly indicates that it would be inappropriate to consider 
the ARCTAS-B data as representative to the average atmo- 
spheric state of Western Arctic in July 2008. 

The regional mean tagged CO vertical distributions shown 
in Fig. 8a for April 2008 indicate that the Western Arctic 
air was largely impacted by ASFF from the surface to the 
upper troposphere in addition to background CO level repre- 
sented by the “Other CO”. This Asian contribution exceeds 
the combined contributions from NAFF and EUFF, partic- 
ularly in free troposphere. The contribution from biomass 
burning sources (BOBB + NBBB) is about 2/3 of the Asian 
fossil fuel contribution in April. Our results support the find- 
ings that Asian anthropogenic emissions are the most impor- 
tant source to the enhanced Arctic CO in the troposphere in 


April (Fisher et al. 2010; Wameke et al., 2010). Neverthe- 
less, the background CO, which accounts for about 1/3 of the 
Western Arctic CO with the majority coming from the CH4 
oxidation, is still the dominant source to the mean CO level 
in the Arctic. The standard deviation (with respect to spatial 
variation) is very small (within 15 ppbv or 10 % of total CO) 
in the lower Arctic atmosphere, indicating a relatively homo- 
geneous CO distribution. The European pollution, which had 
been a major source of Arctic pollution in the earlier decades 
(Quinn et al., 2007, 2008; Shaw 1995), contributes less than 
10 % to CO concentrations in April 2008, and is mostly con- 
fined to the lower and middle troposphere. Contribution of 
NAFF is about half of ASFF throughout the troposphere. 

In July 2008 (Fig. 8b), ASFF and BOBB were the most 
important sources for Arctic pollution (about 35 % and 25 % 
of non-background CO). EUFF contribution was very low 
(less than 5 %). Previous studies indicate that European pol- 
lution is largely enhanced during the high phase of the North 
Atlantic Oscillation (NAO + ) compared to its low phases 
(NAO - ) (Eckhardt et al., 2003). The NAO index was —1.07, 
— 1.73, —1.39, and —1.27 in April, May, June, and July, 
2008, respectively (http://www.cpc.ncep.noaa.gov/products/ 
precip/CWlink/pna/nao.timeseries.gif), suggesting a weaker 
than normal northward transport from Europe to the Arctic. 

The contribution of ASFF to the Arctic mean CO is much 
less in summer than in April, which is due to a lower CO 
concentration and a weaker poleward advection in summer 
than in spring. There is a CO enhancement in the middle to 
upper troposphere around 6-9 km, which suggests that pol- 
lution is transported to the Arctic at high altitudes after being 
lifted over the source regions in the midlatitudes by strong 
convection in July. 

We sort the GEOS-5 total CO according to the dominant 
CO tag for each grid box within the Western Arctic domain, 
as explained in Sect. 3.3. The probability distribution func- 
tion (PDF) of modeled total CO was calculated for each of 
the five source groups, four of which are shown in Fig. 9. 
The EUFF tag was excluded from the figure since the sam- 
pled total CO was almost never dominated by the tagged Eu- 
ropean pollution. The PDFs from the observation and model 
along flight tracks are also shown in the figure as an addi- 
tional model evaluation. 

Figure 9 implies qualitatively how polluted Arctic air 
would be due to each outstanding source and type. During 
ARCTAS-A the Arctic air would be very clean if the major- 
ity of pollution came from NBBB (i.e., the probable total CO 
mixing ratio is very low). In this case the air from the non- 
boreal region must travel long distances before arriving in the 
Arctic. During ARCTAS-B the CO PDFs show a longer tail 
toward higher CO compared to the spring distributions. This 
tail results mostly from the biomass burning sources (BOBB 
and NBBB) with high-intensity episodes (CO > 3 50 ppbv), 
or from a relatively local anthropogenic source (NAFF) with 
frequently high Arctic CO (CO ^200 ppbv). In addition, the 
most probable CO from BOBB during July was less than 


Atmos. Chem. Phys., 13, 4707-4721, 2013 


www.atmos-chem-phys.net/13/4707/2013/ 


H. Bian et al.: Source attributions of pollution to the Western Arctic 


4717 


April July 



0 100 200 300 400 500 0 100 200 300 400 500 



0 100 200 300 400 500 0 100 200 300 400 500 

CO (ppb) CO (ppb) 


Fig. 9. The probability distribution function (PDF) of total CO over 
two ARCTAS campaign periods and over the 4 tagged CO cate- 
gories. Please note EUFF is not shown since CO from European 
anthropogenic emissions seldom dominated the Arctic CO during 
the campaign. The green areas are the PDFs of GEOS-5 CO over 
the Arctic region 50N-90N and 190E-320E and up to 200 hPa. The 
blue and red lines are the PDFs of DC-8 and GEOS-5 sampled along 
flight tracks and N is number of 1 min data point. 


during April although both local and Russian BB CO emis- 
sions during summer were larger than during spring. This 
indicates that a smaller fraction of Russian BB CO was trans- 
ported to the Arctic in summer than in spring. 

A campaign objective was to estimate the emissions from 
boreal forest fires and to investigate the near-field chemical 
evolution of the fire plume. MODIS data over Siberia from 
2000-2009 indicated that fire counts during April 2008 were 
the highest recorded during that period, while the summer 
peak was second only to 2003 (Singh et al., 2010). Canadian 
fires during ARCTAS-B were near their normal level (Soja 
et al., 2008), while California experienced one of the largest 
episodes of wildfires in summer 2008, with over 1000 fires 
(Singh et al., 2010). Episodes from all these biomass burn- 
ing emissions were detected during ARCTAS-B, as shown 
in Fig. 9. In April the Siberian fires were measured over 
Fairbanks, so BOBB gives the highest probable CO. During 
July, California fires (i.e., NBBB) are most prominent over 
Canada, as shown by the longer tail toward high CO. Lo- 
cal Canadian fires (i.e., BOBB) can cause extremely polluted 
episodes as shown by the small bumps at extremely high CO. 


5 Conclusions 

We have analyzed CO and BC measured by the NASA DC -8 
aircraft during the 2008 ARCTAS spring and summer cam- 
paigns. These data were used to evaluate the GEOS-5 model 
simulations. We have also used the observed CH 2 CI 2 , a fossil 
fuel tracer, and CH 3 CN, a biomass burning tracer, to evaluate 
the capability of tagged CO tracers in characterizing air mass 
origins, and the ratio of BC / CO to estimate the age of air. Fi- 
nally, the model is used to attribute the pollutants regionally 
in the Western Arctic to their source origins. 

The model results agree well with the aircraft mea- 
surements along the flight tracks, particularly representing 
long-distance transport of air pollution tracers (from Asian 
sources) in April. The model reveals that the ARCTAS DC -8 
measurements are representative of the regional Arctic pol- 
lution in the spring (April, ARCTAS -A) due to relatively ho- 
mogeneous tracer distribution. The aircraft data alone, how- 
ever, are insufficient to provide a comprehensive and rep- 
resentative picture of Arctic pollution in the summer (July, 
ARCTAS-B) because the flights targeted local fire plumes. 
The tagged CO tracers are able to characterize the air mass 
source regions, as shown by the clear correlation between 
ASFF CO and CH 2 C1 2 and between BOBB CO and CH 3 CN. 
The capability of tagged CO tracers in identifying pollution 
origins is further demonstrated in a case study from 8 July. 

Our model results indicate that on average in both cam- 
paign periods the Western Arctic CO is dominated by back- 
ground CO produced by CH 4 oxidation. Asian anthropogenic 
pollution stands for the largest foreign source for the en- 
hancement of the Western Arctic CO in April (accounting 
for ~25 % of the total Arctic CO), during which biomass 
burning contributes roughly 2/3 as much as ASFF. In July 
both ASFF and BOBB have comparable contributions to the 
mean CO. Biomass burning also makes a large contribution 
to the Arctic pollution variability in July. On the other hand, 
we found that European sources seldom made important con- 
tributions to the CO in the campaign domain and period, per- 
haps because the NAO was negative. 

Our results show that transported pollution from Asian 
emissions is the most important foreign fossil fuel and 
biomass burning sources of CO and BC to the Western Arctic 
throughout the troposphere during April, and in the middle- 
upper troposphere (6-10 km) during July, with stronger 
transport capability during spring than summer. 

The ratio of BC / CO can be used as an indicator of air 
mass age and origin. The air masses originating from Asian 
anthropogenic emissions were relatively old, while those 
from boreal biomass burning emission were relatively young 
in both campaigns. The ratio varies most in the July BOBB 
case due to the two distinct origins from local Arctic boreal 
fires and from Eurasia. The ratio of ASFF is larger in April 
than in July, which is consistent with more BC wet scaveng- 
ing occurring during the Asia-Arctic transport in July. 
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During ARCTAS-B the local boreal biomass burning 
sources from the standard QFED seemed to significantly 
overestimate the emissions. The ARCTAS observations 
are used to constrain this emission. Even so, the model- 
measurement comparisons indicate that the model still has 
difficulties to capture the location and time of local boreal 
forest fire emissions. Despite this deficiency the model is use- 
ful in representing the Western Arctic region-wide pollution 
as the regional mean CO is less impacted by the local events. 
Furthermore, using a varying conversion time for BC from 
hydrophobic to hydrophilic (i.e., 2.5 days in ARCTAS -A and 

1 .25 days in ARCTAS-B) yields a better model-observation 
agreement of BC mass concentration. 

Supplementary material related to this article is 
available online at: http://www.atmos-chem-phys.net/13/ 
4707/2013/acp-13-4707-2013-supplement.pdf. 
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